function delta_min=findDelta(sites, tour)

[val, Isite, Itour]=intersect(sites(:,1:2), tour(:,1:2), 'rows');
delta_z=tour(Itour,3)-sites(Isite,3);
delta_min=min(delta_z);